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Probes that measure the local thermal properties of systems out of equilibrium are emerging as 
new tools in the study of nanoscale systems. One can then measure the temperature of a probe that 
is weakly coupled to a bias-driven system. By tuning the probe temperature so that the expecta¬ 
tion value of some observable of the system is minimally perturbed, one obtains a parameter that 
measures its degree of local statistical excitation, and hence its local heating. However, one antici¬ 
pates that different observables may lead to different temperatures and thus different local heating 
expectations. We propose an experimentally realizable protocol to measure such local temperatures 
and apply it to bias-driven quantum dots. By means of a highly accurate open quantum system 
approach, we show theoretically that the measured temperature is quite insensitive both to the 
choice of observable and to the probe-system coupling. In particular, even with observables that are 
distinct both physically and in their degree of locality, such as the local magnetic susceptibility of the 
quantum dot and the global spin-polarized current measured at the leads, the resulting local temper¬ 
atures are quantitatively similar for quantum dots ranging from noninteracting to Kondo-correlated 
regimes, and are close to those obtained with the traditional “local equilibrium” definition. 

PACS numbers: 72.20.Pa, 79.10.-n, 71.27.-l-a, 73.63.Kv 


Temperature is a thermodynamic quantity of funda¬ 
mental importance in the description of systems at equi¬ 
librium. However, the extension of this concept to sys¬ 
tems far from equilibrium is not obvious. From an op¬ 
erational point of view, temperature is defined as the 
quantity measured by a thermometer coupled to the sys¬ 
tem with which it reaches thermal equilibrium. Precisely 
because the thermometer plus the system reach a state 
of global equilibrium when coupled, the quantity that is 
measured by the thermometer is then attributed to the 
system in the limit of weak coupling and negligible heat 
capacity of the thermometer. This definition is no longer 
valid for a system out of equilibrium such as one driven 
by a constant bias [1]. In this case, electron-electron in¬ 
teractions and electron-phonon interactions are expected 
to induce electronic and ionic “temperatures” different 
from those of the same system at equilibrium [1]. The 
question then is: What are these temperatures and how 
do we measure them directly! 

Several solutions have been proposed. For instance, 
Engquist and Anderson have introduced the concept of 
ideal potentiometer/thermometer [2]. In this case, local 
chemical potential and local temperature are defined by 
a “local equilibrium condition”: the net particle and heat 
current flowing through the potentiometer/thermometer 
are set to zero [3-11]. Although such a definition is ap¬ 
pealingly intuitive and has been used extensively in the 
past, we do not have means to directly measure heat cur¬ 
rents (unlike particle currents for which ammeters are 
available) [12], and therefore its experimental realiza¬ 


tion is very limited. Other theoretical definitions include 
the use of approximate distribution functions [13], infor¬ 
mation compressibility [14], or generalized fluctuation- 
dissipation theorems [15], to name just a few (for a more 
complete set of definitions see, e.g.. Ref. 12). All of these, 
however, suffer from some limitations in their experimen¬ 
tal realizations and therefore are also of limited use. 

On the experimental side, methods have been devised 
to determine local temperatures by monitoring proper¬ 
ties that are sensitive to thermal fluctuations, such as 
the bond rupture force [16], the junction lifetime [17, 18], 
the mechanical stretching distance [19], the bias-driven 
current noise [20], and the surface-enhanced Raman in¬ 
tensities [21, 22]. These probes, however, provide only 
an indirect measurement of a “local temperature” leav¬ 
ing the original question still open. 

There is, instead, an increasing body of experimental 
studies in which thermal probes are coupled directly to 
driven nanoscale systems with a resolution that is ap¬ 
proaching hundreds of nanometers or less [23-30], thus 
making them ideal as local thermometers. With these 
thermometers an operational definition of temperature 
has been proposed by imposing a minimal perturbation 
condition [31], in which the temperature of the probe, 
Tp, is varied, while monitoring some observable of the 
system, in such a way that the expectation value of that 
observable is minimally perturbed. The probe temper¬ 
ature satisfying this condition is then a parameter at¬ 
tributed to the system, which characterizes its local ex¬ 
citations out of equilibrium. This type of definition has 
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been used, for instance, in the study of the thermoelectric 
response of nanoscale systems to applied thermal gradi¬ 
ents, leading to the prediction of temperature oscillations 
[10, 32], It is an operational definition that is relatively 
straightforward to implement experimentally. However, 
it leaves open the prescription of what type of observable 
one should use, and whether different observables lead 
to quantitatively different temperatures. In addition, in 
the original paper [31], the thermal probe was assumed 
to be a bosonic system, and hence no electric current 
could flow between the system and the probe. This is 
quite a strong limitation since the experimental probes 
that are being developed are generally fermionic systems 
[23-30]. In this case, an extra constraint related to the 
local chemical potential of the thermal probe needs to be 
introduced. 

In this paper, we discuss a protocol to measure the 
local temperature of a system out of equilibrium that 
is coupled to a thermal probe as schematically shown in 
Fig. 1. We consider a fermionic probe and allow the mea¬ 
surement of an observable of the system in such a way 
that the minimal perturbation condition be satisfied. We 
show that by appropriately setting the chemical poten¬ 
tial of the probe at the beginning of the measurement, 
the choice of observables quite distinct both physically 
and in their degree of locality lead to local temperatures 
that are quantitatively similar even for strongly corre¬ 
lated systems, and very close to those obtained with the 
often-used local equilibrium temperature. These results 
lend support to the extrapolated parameter as the “tem¬ 
perature” of the system. 

To be specific, we consider a quantum dot (QD) in 
contact with two leads as sketched in Fig. I. Under a bias 
voltage or a thermal gradient across the two leads, the 
temperatures (chemical potentials) of left and right leads 
are T]^ and Tfj (/i^ and iir), respectively. To determine 
the local temperature T* and local chemical potential 
of the QD, the dot is coupled to a third lead (the 
probe), whose chemical potential jip and temperature Tp 
are tunable. 

The first step of the protocol we suggest is the determi¬ 
nation of the chemical potential of the probe. Ideally, at 
zero bias (equilibrium) the minimal perturbation condi¬ 
tion [31] should yield exactly the background equilibrium 
temperature Teq. This can be accomplished as follows. 
In the presence of an applied bias voltage or a thermal 
gradient [33], we first determine /r* as 

( 1 ) 

The weight coefficients Cl and Cfl are determined by [34] 

^ _ 1 _ _ IpjTg, p-Q.) _ 

Ip{n,^iR)-Ip{TR,^^R) ■ 

Here, Ip{Tg,fig) is the electric current measured at the 
probe, by setting the chemical potential and temperature 
of the probe to be identical with those of lead-a. 


detector 



FIG. 1. Schematic diagram of the proposed protocol for the 
measurement of local temperature in a QD connected to two 
leads (L and R). A weakly coupled probe with tunable fip 
and Tp is used. The expectation value of an observable. Op, 
is monitored while Tp is varied. See the main text for details. 

We then set fj,p = , and monitor the change of a 

given system observable O = (O) as Tp varies. The lo¬ 
cal temperature T* is finally determined by the minimal 
perturbation condition [31] 

T* = argmin |C>o -b 50p - Op{Tp, n*)\. (3) 

Tp 

Here, Oq is the expectation value of O measured without 
the probe, while Op[Tp,^*) is its measured value with 
the probe coupled to the dot. The nonzero probe-dot 
coupling, Ap, results in a finite perturbation to the in¬ 
trinsic dot properties. This effect is accounted for by the 
correction term 50p in Eq. (3) determined by 

^Op = Cl Op{Tl, /tl) -b Op{TR, fi^) — Oq. (4) 

Here, Op{Ta, fig) is measured by setting Tp and /ip to be 
identical to their counterparts of lead-a. 

Note that the protocol proposed above is easily realiz¬ 
able experimentally, and is universally applicable to any 
QD system. Also, for all the QDs studied in this work, 
the minimal perturbation of (O) searched for in Eq. (3) is 
always found to be zero, i.e., Op(Tp, /i*) = Oq + SOp can 
always be satisfied at a certain Tp. This indicates that 
the chosen observable Op changes monotonically and sen¬ 
sitively as Tp varies in the vicinity of T*. 

It is important to verify that at exactly zero bias the 
local temperature measured by the minimal perturbation 
condition is exactly the physical equilibrium temperature 
Teq. This can be seen from Eqs. (I)-(4). At zero bias, 
we find that from Eq. (I) the local chemical potential is 
just the equilibrium Eermi energy. Then Eq. (3) reduces 
to the simple form Op(Tp) = Op(Teq), with the trivial 
solution Tp = Teq. 

We now apply Eqs. (I)-(4) to QD systems described by 
a single-impurity Anderson model (SIAM) ]35, 36]. The 
total Hamiltonian \s H = FIdot + FFiead + Flcoupiing- The 
dot is described by iFdot = idfid + Here, fid = 
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where aj (os) creates (annihilates) a 
spin-s electron on the dot energy level e^, and U is the 
on-dot electron-electron interaction strength. TJiead = 

^2aks ^aks ^aks and i?coupling = X/afes ^aks + 

H.c. represent the noninteracting leads and dot-lead cou¬ 
plings, respectively. Here, (daks) creates (annihi¬ 

lates) an electron on the orbital \k) of lead-a {a = L,R 
or p); and tak is the coupling strength between the dot 
level and lead orbital |fc). 

We employ a hierarchical equations of motion (HEOM) 
approach to compute the reduced density matrix (p) 
of open fermionic systems [37, 38], so as to character¬ 
ize the equilibrium and nonequilibrium properties of the 
SIAM [39, 40]. The effect of the leads is captured by the 
hybridization functions Aq,(w) = {tak^^die — Cak)- 

For numerical convenience, a Lorentzian form AQ(a;) = 
/[{uj — PaY + W‘^] is adopted, where A^ is the cou¬ 
pling strength between the dot and lead-a, and W is the 
lead band width. Hereafter, A = A^-I-Afl is taken as the 
unit of energy. The HEOM approach has been used to 
study static and dynamic Kondo effects in QDs [41-45], 
and it is in principle exact if all orders of hierarchy ex¬ 
pansion were included. In practice, the numerical results 
converge to the exact values rapidly and uniformly with 
the increasing truncation level of hierarchy Atrun- Once 
the convergence is achieved, the results are guaranteed 
to be quantitatively accurate [38]. 

We compute the system observable via O = tr(pO), 
as well as the energy distribution of electric and heat 
currents flowing into lead-a as required to determine T* 
with the “local equilibrium condition” [46-49] 

X + 2ifp(uj - Pa)lm[0’'(a;)]} . (5) 

Here, we have set e = h = 1; k = 0 and 1 correspond 
to the electric and heat currents, respectively; fpioj) is 
the Fermi function; and the lesser and retarded Green’s 
functions and are computed from correlation func¬ 
tions. The global electric and heat currents flowing into 
lead-a are da = J dojj^{uj) and Ja = J doo respec¬ 
tively. Note that the electron-phonon interactions and 
the phonon contribution to heat current have been ne¬ 
glected, since their effects are negligibly small in QDs at 
low temperatures [50]. 

We consider first a noninteracting QD {U = 0) under 
a bias voltage pa = —pL = 5 ^- The HEOM method 
has the virtue that, for noninteracting systems the hi¬ 
erarchy terminates automatically at Atrun = 2 without 
any approximation [37]. Figure 2(a) plots T* deter¬ 
mined by Eqs. (l)-(4) versus the background tempera¬ 
ture T = Tl = Tr. We have set the dot-lead couplings 
to be asymmetric with A^, = The weight coeffi¬ 

cients in Eq. (1) are then C,a = Aq,/A [34], and the local 
chemical potential is thus p* = ^ pL = 



FIG. 2. (a) Local temperature T* of a noninteracting QD 
determined by different protocols versus background T. The 
dashed line is a guide to eyes.(b) Variation of the local mag¬ 
netic susceptibility and total spin current /)(* with Tp 
at a background T — 0.133A. The parameters adopted are 
(in units of A): td = —2.67, [7 = 0, Al = 2A_r = 0.67, 
PR = —pL = 0.53, and W = 40. 


To have an accurate measurement of T* and thus 
test its robustness with respect to the choice of observ¬ 
ables, the probed observable O must vary sensitively with 
Tp. We then choose some spin-related properties, be¬ 
cause spin polarization processes require low excitation 
energies. These include the local magnetic susceptibility 
X™ = * fo d7([ffiz(f),ffiz(0)]) with = gpB(nt - n;)/2, 
and the global spin-polarized current 7™ = (/q^) — (lai) 
with las = i ^coupling]- Here, g is the gyromag- 

netic ratio and ps is the Bohr magneton. The different 
degree of locality of these quantities provides an even 
stronger test for our protocol. As shown in Fig. 2(b), 
both y™ and I™ vary monotonically with Tp, and their 
perturbations by the probe reduce to zero at roughly 
the same Tp ~ 0.28A. Therefore, from Eq. (3) we have 
T* = 0.28A at T = 0.133A. Figure 2(a) compares T* 
obtained by different protocols over a large range of T. 
In the same plot we also show that T* determined by the 
minimal (zero) perturbation of y™ or J™ agrees closely 
and consistently with that obtained with the local equi¬ 
librium condition. A closer look at both the electric and 
heat currents that flow between the probe and the system 
when the minimal perturbation is satisfied reveals that 
they are both close to zero [51], thus explaining the agree¬ 
ment between the two protocols. It is also found that, 
under a bias voltage, T* is always higher than T, and 
their deviation increases as T decreases. In particular, 
T* maintains a finite value (0.25A) even as T —0. This 
indicates that the local heating feature becomes more 
visible at a lower background T [1]. 

We now investigate interacting QDs under an antisym¬ 
metric bias voltage. We choose to examine a half-filling 
dot with U = —2ed = 2.4A. It has been shown (through 
the temperature-dependent conductance) that this QD 
exhibits prominent Kondo features at T < Tr- [38], with 
Tk = 0.82A being the characteristic Kondo temperature 
[52]. This is evident from the inset of Fig. 3(a) where 
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FIG. 3. (a) Local temperature T* of an interacting QD de¬ 
termined by different protocols versus background T. The 
dashed line is a guide to eyes. The inset shows the Kondo 
spectral peak of the dot at various T. (b) Variation of system 
observables x"* a-nd with Tp, at a background T = O.IA 
lower than the Kondo temperature Tk = 0.82A. The pa¬ 
rameters adopted are (in units of A): Sd = —1.2, U = 2.4, 
Al = An = 0.5, fj,n = -^.l = 0.25, and W = 5. 
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FIG. 4. Local temperature T* of a noninteracting QD de¬ 
termined by different protocols versus the system-probe cou¬ 
pling strength Ap. The background temperature is fixed at 
T = 0.133A. The parameters adopted are (in units of A): 
f-d = —2.67, U = 0, An = 2Aij = 0.67, jj.R = —fJ-n = 0.53, 
and W = 40. 
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the dot spectral function A{u}) = — ;hlin[5’'(a;)] is plotted 
for various background T showing that the QD exhibits 
more prominent Kondo feature as T lowers, as confirmed 
by the higher and sharper Kondo spectral peak centered 
at w = /r*. The chemical potentials of the leads are 
fiR = —fin = 0.25A. With the symmetric dot-lead cou¬ 
plings {An = An), the local chemical potential on the 
dot is n* =0 by Eq. (1). 

Figure 3(a) compares T* determined by different pro¬ 
tocols. Limited by computational resources at our dis¬ 
posal, the HEOM calculations are done at iVtrun = 4 for 
the interacting QD [53]. The results are considered to 
be highly accurate within the explored range of temper¬ 
atures, since the computed conductance well reproduces 
the Kondo scaling relation [38, 45]. As shown in Fig. 3(a), 
the local temperatures T* determined by the minimal 
(zero) perturbation of y™ and agree remarkably well 
with each other, as well as with the local equilibrium 
temperature. This is again because the electric and heat 
currents flowing into the probe are both close to zero 
when the minimal perturbation is satisfied [51]. 

Similar to the noninteracting QD, the local tempera¬ 
ture T* on the interacting dot is always higher than the 
background T, and it approaches a finite value (T* ~ 
0.15 A) as T —7> 0. As clearly indicated by Fig. 3(b), while 
both and 7™ vary monotonically with Tp, the lat¬ 
ter has a much more sensitive temperature dependence. 
More importantly, the perturbations of both y™ and 7£* 
reduce to zero at almost the same Tp, which highlights 
again the generality of our operational definition. 

Finally, we investigate the influence of the coupling 
strength between the probe and the system. For numeri¬ 
cal convenience we choose to examine the noninteracting 
QD explored in Fig. 2. Figure 4 shows T* determined by 
using the minimal perturbation as well as the local equi¬ 
librium conditions. As the probe-dot coupling reduces to 
zero, all these temperatures approach a constant value as 


expected. More telling though is the fact that, even with 
a relatively large probe-dot coupling strength (compara¬ 
ble with the dot-lead couplings) the resulting T* have 
only a minor change of ^ 7%, indicating again the ro¬ 
bustness of the proposed protocol. We thus conclude that 
the measured local temperature depends rather insensi¬ 
tively on the probe-dot coupling Ap, which is favorable 
for experimental realizations. 

In summary, we have proposed an operational def¬ 
inition of local temperature for bias-driven QDs us¬ 
ing a “minimal perturbation condition”, as represented 
by Eqs. (l)-(4). The same definition is also applica¬ 
ble to QDs subjected to external thermal gradients [54]. 
The operational definition applies equally well to sys¬ 
tems ranging from noninteracting to Kondo-correlated 
regimes. Since this definition does not require measure¬ 
ments of heat currents, its experimental realization is 
straightforward. The “minimal perturbation condition” 
thus provides a useful practical means to examine local 
electron excitations in a nonequilibrium process, in which 
local heating plays an important role. 
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I. REMARKS ON THE WEIGHT COEFFICIENTS Cl AND Ci^ 


Suppose Tp = Tr. The electric current between left and right leads is driven only by the applied 
bias voltage (difference in chemical potentials). With the chemical potential of probe Hp aligned 
with the dot local chemical potential pp, the net electric current flowing into the probe is exactly 
zero. This is because the electric currents coming from the left and right leads cancel out exactly: 


{p* - pl)Glp = {pR - p*)GRp. 


(SI) 


Here, Gap (a = L, R) is the electron conductance between lead-a and the probe. Rearranging 
Eq. (SI), we have 


T = 


Grp 




PR + 


Glp 




-Ur.. 


(S2) 


By setting pp = pa and Tp = Ta, the probe can be deemed as a part of lead-a. Therefore, the 
electric current flowing through the probe Ip amounts to 


IpiTp, pl) = Grp {pL - Pr), 

Ip{TR, pr) = Gpp {pR - Pl), (S3) 


G I (T LL 'l 

and hence comparing Eq. (S2) and Eq. (1) of main text, it is clear that the 

weight coefficient (a should be proportional to the two-terminal conductance Gap, as follows 


Ca = 


G, 


ap 


Glp + Grp 


= 1 - 


IpiTa, Pa) 


Ip{TL, Pl) - Ip{TR, Pr) 


(S4) 


Here, we have assumed that Gap does not change significantly with pp and Tp. Obviously, the 
weight coefficients should normalize to unity. Cl + Cr = 

Eor the single-impurity Anderson model (SIAM) investigated in the main text, if all leads have 



2 


the same form of hybridization function, i.e., Aa{uj) = Aar]{Ld) {a = L, R and p), we have 


Ip{TL,fJ-L) = ^ 


A, 


Al +A, 


2\ (Al + Ap)A/{ 


TT J Al + Aji + Ap 


X du - ^l) - - fiR)] Im[^^(a;)], 


(S5) 


Ip(TR, pr) = 


A_R + Ar 


{Ar + Ap)Al 

TT J Ai + A^ + Ap 


X du; [ffinioj - Pr) - - pl)] Im[a’'(w)], 


(S6) 


with f3a = l/7h. Here, the probe is deemed as part of lead-a, and the SIAM amounts to a two- 
terminal model. Therefore, Eq. (9) of Ref. 1 can be used. Equations (S5) and (S6) immediately 
lead to the relation of 

Ip(Tl,Pl) ^ 

Ip{Tr,pr) Ar- ^ > 

Equation (S7) thus offers a practical means of determining the ratio Ar/Ar in experiments. Con¬ 
sequently, by combining Eqs. (S4) and (S7), we have 


Ca 


A 


a 


Al -I- Ar 


(S8) 


for a = L or R. 

Since the electric current flowing into the probe Ip can be measured straightforwardly, the 
weight coefficients {Ca} are obtained readily from Eq. (S4). Table SI lists the weight coefficient (r 
for the quantum dots (QDs) studied in Figs. 2 and 3 in the main text. The numbers obtained via 
Eq. (S4) and Eq. (S8) agree precisely with each other. 



Fig. 2 

Fig. 3 

Cl 

noninteracting QD 

Kondo-correlated QD 

Eq. (S4) 

0.6666 

0.5000 

Eq. (S8) 

2 

3 

1 

2 


Table SI: The weight coefficient (r for the QDs studied in Figs. 2 and 3 in the main text. The numbers are 
obtained via Eqs. (S4) and (S8), respectively. 


As demonstrated clearly in the main text, the local temperature measured by using the minimal 
perturbation condition agrees closely and consistently with that determined by the local equilibrium 
condition. This is because both the electric and heat currents flowing into the probe are close to 
zero when the minimal perturbation condition is satisfied. Here, we verify that the values of both 
Ip and are indeed negligibly small when the minimal perturbation is reached. The following 
Tables S2 and S3 concern the noninteracting and interacting QDs investigated in the main text, 
respectively. It is clearly seen from the tables that, even after being scaled by the small dot-probe 
coupling strength Ap, the electric and heat currents flowing through the probe are still much 
smaller than the currents flowing into the left and right leads. 
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scaled electric current (e/h) 

scaled heat current (A/h) 

T/A 

Ip/Ap 

Il/Al 

Ir/I^r 


^?/Al 

JI/Ar (JH + JH)/A 

0.133 

2.9 X 10-3 

-0.37 

0.74 

-9.8 X 10- 

-3 0.16 

0.47 0.26 

0.667 

-5.5 X 10-2 

-0.54 

1.08 

-0.17 

-0.25 

1.65 0.39 


Table S2: Electric and heat currents flowing into leads coupled to a noninteracting QD when the minimal 
(zero) perturbation condition is satisfied. The parameters adopted are (in units of A): = —2.67, U = 0, 

Al = 2A/f = lOOAp = 0.67, fXR = —/xl = 0.53, and W = 40. 



scaled electric current (e/h) 

scaled heat current (A/h) 

T/A 

Ip/Ap Il/Al 

Ir/^r 

J^/Ap JI/Al 

JI/Ar (Jl + Jl)/A 

0.4 

-1.3 X 10-® -1.32 

1.32 

-6.9 X 10-3 1.33 

1.33 2.66 

0.8 

-4.9x10-33 -1.01 

1.01 

-3.9 X 10-3 1.02 

1.02 2.04 


Table S3: Electric and heat currents flowing into leads coupled to an interacting QD when the minimal 
(zero) perturbation condition is satisfied. The parameters adopted are (in units of A): ej, = —1.2, U = 2.4, 
Al = Ar = lOOAp = 0.5, fiR = -^iL = 0.25, and W = 5. 

II. APPLICATION OF LOCAL EQUILIBRIUM CONDITION 

Thermodynamically, the dot can be deemed as a reversed heat engine, and the leads act as heat 
baths as well as electron reservoirs. The internal energy of the dot is 

U = Q + W + iJ.*{nd), (S9) 

with Q and W being the heat and work gained by the dot, respectively. Taking the time derivatives 
of both sides of Eq. (S9) leads to 

J| = + P + (SIO) 

Here, Jj" is the energy (heat) current flowing into the dot, and P is the electric power. In a 
stationary state, U is a constant, and thus = d{hd)/dt = 0. Therefore, P = —J^ = Yla 

Figure SI (a) and (b) depict the energy distribution of bias driven electric and heat currents 
and Jq(w)] in the noninteracting QD system studied in Fig. 2 of main text, respectively. All 
data are computed with the hierarchical equations of motion (HEOM) approach. In the absence 
of the probe, j^{oj) = holds for any cu, which ensures the conservation of electric current 

= 0). In contrast, the heat currents flowing through left and right leads do not cancel out 
(^Q, / 0) - the dot behaves as a hot spot and dissipates heat into both leads. 

To utilize the local equilibrium condition, the chemical potential and temperature Tp of 
the probe are tuned until the electric and heat currents flowing into the probe [Ip = J dui jp{uj) 
and = f dujjp(uj)] are both zero. Figure SI (c) and (d) display jp(io) and jp(Lo) when local 
equilibrium is reached (Ip = = 0), respectively. 
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co/A ca/A 



co/A 


co/A 


Figure SI: Panels (a) and (b) depict the energy distribution of electric and heat currents flowing into a-lead 
(a = L,R), respectively. Panels (c) and (d) plot the energy distribution of electric and heat currents flowing 
into the probe when the local equilibrium condition Jp = = 0 is satisfied, respectively. The parameters 

adopted are (in units of A): ed = —2.67, U — 0, = 2Aii = lOOAp = 0.67, = 0.53, W = 40, 

and Tr = Tr = 0.133. 


III. LOCAL TEMPERATURE OF A QUANTUM DOT SUBJECTED TO A THERMAL 

GRADIENT 

We now demonstrate that the proposed operational definition of local temperature, Eqs. (1)- 
(4) of main text, is also applicable to QDs subjected to an applied thermal gradient, i.e., when 
AT = Tr — Tr ^ 0. Note that for this to be true, the bias voltage must be zero {fiR = hr), and 
thus the electric current is driven only by the thermal gradient AT. 

In the absence of bias voltage, Eq. (1) of main text becomes trivial since fj,* = fiR = fiR. 
Therefore, the chemical potential of probe is fixed at Hp = m = fJ-R, and we only need to tune Tp. 
When Tp = T*, the electric current flowing into the probe vanishes. We thus have 

{T*-Tl)lI^ = {Tr-T*)LI^. (Sll) 

Here, L^p is the thermoelectric transmission coefficient between lead-a and probe. Through an 
analysis similar to that in Sec. I, one arrives at a simple expression of T* as 

T* = ClTl + (:rTr. (S12) 

Here, the weight coefficients Cl and Qr are also given by the last equality of Eq. (S4). However, 
Eq. (S4) turns out to be a rather crude estimate of T*, as shown in Eig. S2 below. The accurate 
measurement of T* is achieved by using Eqs. (3)-(4) of the main text. 
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For instance, consider a general case that dot-lead couplings are asymmetric with = 5Aji. 
The weight coefficients Cl = § aiid Cfl = g agree well with the “measured” values of Cl = 0.8333 
and Cr = 0.1667 by Eq. (S4). Figure S2 compares T* obtained by different approaches. Clearly, T* 
determined by minimal (zero) perturbation of x”* and I™ agree remarkably well with the prediction 
of local equilibrium condition. 



Figure S2: T* of a noninteracting QD versus the temperature of left lead T^. The parameters adopted are 
(in units of A): = —3.33, U = 0, = SA^^ = 0.83, = = 0, AT = Tr — Tr = 0.17, and W = 50. 
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